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Abstract 

The  Wireless  Integrated  Microsystems  (WIMS)-environmental  monitor  testbed  (EMT)  is  a  multi-component  microelectromechanical 
system  (MEMS),  incorporating  complementary  metal  oxide  semiconductor  (CMOS)  materials  for  high-precision  circuits  used  for  integrated 
sensors  such  as  micro-g  accelerometers,  micro-gyroscopes,  and  pressure  sensors.  The  WIMS-EMT  duty  cycle,  like  many  autonomous  MEMS 
systems,  has  low-power  standby  periods  for  sensing,  and  high-power  pulses  for  R/F  transmission  and  reception.  In  this  paper,  we  present 
results  of  three  strategies  for  providing  power  to  this  system,  including  (1)  specification  of  a  single,  aggregate  power  supply,  resulting  in  a 
single  battery  electrochemistry  and  cell  size;  (2)  specification  of  several  power  supplies,  by  a  priori  division  of  power  sources  by  power  range; 
and  (3)  specification  of  an  arbitrary  number  of  power  “bundles,”  based  on  available  space  in  the  device.  The  second  approach  provided  the 
best  results  of  mass  (0.032  kg)  and  volume  (0.028  L)  among  the  three  approaches.  The  second  and  third  approaches  provided  the  best  battery 
lifetime  results;  both  systems  produced  lifetimes  in  excess  of  2E3  h.  Future  work  will  incorporate  CMOS  operational  amplifier  (op-amp) 
technologies  to  accommodate  large  voltage  fluxes  in  many  MEMS  devices,  and  implementation  of  our  approaches  into  a  user-friendly  code. 
©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  multicomponent  Wireless  Integrated  Microsystems 
environmental  monitor  testbed  (WIMS-EMT)  is  a  micro¬ 
electromechanical  system  (MEMS)  whose  goal  is  the  re¬ 
alization  of  a  wristwatch-sized  device  capable  of  sampling 
ambient  pressure,  temperature,  humidity,  and  air  quality  [1]. 
The  WIMS-EMT  incorporates  complementary  metal  oxide 
semiconductor  (CMOS)  materials  for  high-precision  circuits 
used  for  integrated  sensors  such  as  micro-g  accelerometers, 
micro-gyroscopes,  and  pressure  sensors. 

Supplying  power  to  CMOS-based  platforms  has  become 
increasingly  challenging,  due  to  nonlinear  scaling  of  power 
consumption  with  respect  to  MEMS  device  size  and  mass, 
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and  difficulties  associated  with  manufacture  of  thin-film  bat¬ 
teries  using  CMOS  fabrication  techniques.  CMOS  devices 
are  also  prone  to  heat-induced  failure  [2-5].  Despite  many 
advances  in  reduction  in  power  consumption  [2,6-8]  of 
novel,  higher-power  and  capacity  micro  and  thin-film  bat¬ 
teries  [9-13]  a  cohesive  methodology  for  design  of  power 
systems  for  MEMS  has  not,  to  our  knowledge,  been  reported. 
Here,  we  introduce  such  an  algorithm,  which  involves  use  of 
system  constraints  on  operating  temperature,  storage  battery 
life,  energy/power  density,  specific  energy/power;  and  bat¬ 
tery  characteristics,  e.g.  rechargeability,  mass,  volume,  and 
lifetime.  CMOS  power  consumption  is  a  function  of  static 
and  dynamic  energy  consumption,  i.e. 

Aotal  =  ^SW  +  PSC  +  P static  +  P\q  ak-  (1) 

Dynamic  switching  power,  Esw,  is  required  to  change 
charge  to  discharge  functions.  Short  circuit  power  consump- 
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tion,  Psc>  results  from  current  traveling  between  the  power 
supply  and  ground  during  gate  switching.  The  static  power 
component,  /Static*  results  from  other  factors,  including  bias 
current.  Finally,  power  consumption  due  to  leakage  cur¬ 
rent,  Pieak  results  from  gate-oxide  tunneling  current,  source 
to  drain  sub- threshold  conduction,  and  pn-j unction  leakage 
[14].  Of  these,  the  dominant  component  of  power  consump¬ 
tion  in  CMOS  devices  is  dynamic  switching  power  (85-90%) 
[2],  and  a  first  order  approximation  of  this  term  can  be  ex¬ 
pressed  as: 

^SW  =  ax Supply/  (2) 

where  x  is  capacitance  due  to  charge/discharge  switching,  a 
is  an  activity-weighting  factor  that  represents  the  probability 
that  a  transition  occurs,  VSupply  is  the  supply  voltage,  and/is 
the  frequency  of  operation.  Many  workers  have  focused  on 
reducing  the  factors  that  most  influence  dynamic  power  con¬ 
sumption,  namely  capacitive  load,  supply  voltage,  switching 
frequency,  and  activity  (e.g.  [2,6-8]).  However,  alteration  of 
some  of  these  parameters  can  independently  reduce  perfor¬ 
mance;  for  example,  reduction  in  supply  voltage  generally 
results  in  slower  speed  of  execution  of  instruction  sets  a  pro¬ 
cessor  can  execute,  bandwidth,  and  clock  speed. 

Reductions  in  dynamic  power  have  been  forecasted,  with 
reductions  in  scale  of  CMOS  technology  to  deep-submicron 
levels.  MOS  transistor  channel  lengths  have  been  produced 
with  dimensions  under  100  nm,  and  device  gate-oxide  iso¬ 
lation  have  been  reported  of  less  than  2  nm.  Thus,  power 
due  to  component  leakage  is  predicted  to  ultimately  surpass 
dynamic  power  [2].  The  ITRS-2001  (International  Technol¬ 
ogy  Roadmap  for  Semiconductors  in  2001),  has  predicted 
that  operating  voltage  will  decrease  from  1.2  V  (2001)  to 
0.5  V  (2013),  and  that  the  maximum  number  of  transistors 
per  chip  of  most  silicon  technologies  will  increase  from  7 
million  (2001)  to  8  billion  (2013)  due  to  circuit  device  tech¬ 
niques  such  as  multi-threshold  CMOS  [15],  adaptive  reverse 
body  biasing  [16],  and  microarchitectural  and  compiler  tech¬ 
niques.  For  example,  wireless  personal  area  communication 
devices  (CD),  high-end  user  interface  devices  (also  called 
master  devices,  MD),  and  low-end  user  devices  are  pro¬ 
jected  to  require  active  power  supplies  in  the  ranges  of  100, 
10,  and  1  mW,  respectively,  and  standby  power  consump¬ 
tion  values  of  10  mW,  1  mW,  and  100  jxW,  respectively  [2,3]. 
Switching  power  consumption  in  these  devices  (using  val¬ 
ues  from  ITRS-2001 — International  Technology  Roadmap 
for  Semiconductors  in  2001),  are  predicted  to  decrease  from 
^300  mW  in  2004,  to  less  than  1  mW  by  2013;  while  static 
power  is  predicted  to  increase  from  0.25  W  in  2004  to  4.25  W 
by  2013  [2,3].  Scaling  of  large  devices  via  implementation 
of  CMOS  technology  will  remain  nonlinear  with  respect  to 
reduction  in  power  supply. 

Specific  energy  and  power  (with  typical  dimensions 
Whkg-1  and  Wkg-1,  respectively)  and  energy  and  power 
densities  (with  typical  dimensions  WhL-1  and  WL-1,  re¬ 
spectively)  are  generally  used  to  characterize  the  mass  and 


volumetric  efficiency  of  a  power  supply  [17],  though  these 
values  alone  are  insufficient  to  prescribe  a  system.  Because 
many  secondary  batteries  exhibit  capacity  fade  and  rate- 
dependence,  along  with  temperature  and  other  environmental 
constraints,  more  information  is  needed  to  select  an  appro¬ 
priate  power  supply. 

Thin-film  electrochemistries  have  been  developed,  of  sizes 
ranging  from  1  to  21 .7  cm2  in  surface  area,  8  [xm  and  3.6  mm 
in  thickness,  with  current  discharge  rates,  voltages,  and 
specific  energy  ranging  from  0.1  to  1.6  C,  1.5-4.0V,  and 
26.3-350  (Whkg-1),  respectively,  listed  in  Table  1  (per  re¬ 
ported  values  in  [9,10,12,13]).  However,  implementation  of 
these  batteries  into  “real-world”  devices  has  been  problem¬ 
atic,  since  CMOS  fabrication  temperatures  routinely  exceed 
250  °C,  which  is  above  the  boiling  point  of  some  battery  ma¬ 
terials  (e.g.  lithium  metal,  at  180.5  °C)  [18].  Tin  oxide  can 
be  processed  and  cycled  at  high  temperatures  [18].  However, 
it  cannot  provide  the  same  capacity  as  lithium  metal,  and  is 
thus  often  an  unsuitable  replacement. 

Factors  influencing  selection  of  power  supply  can  be  sum¬ 
marized  as  follows: 

( 1)  electrochemistry :  cell  potential,  discharge/charge  profile, 
capacity,  and  lifetime 

Conductivity  of  active  materials  largely  determines  the 
maximum  cell  discharge/charge  rates,  and  cell  capac¬ 
ity.  For  example,  Bates  et  al.  [9]  compared  the  perfor¬ 
mance  of  thin-film  lithium  metal  rechargeable  batteries 
which  used  an  amorphous  inorganic  electrolyte,  amor¬ 
phous  V2O5  and  LixM^CU  as  cathode  materials.  They 
found  that  the  performance  of  rechargeable  thin-film 
lithium  batteries  depended  critically  on  methods  used 
to  deposit  the  cathodes,  their  geometry,  operating  tem¬ 
perature,  and  current  density. 

(2)  geometry :  surface  area,  volume,  and  mass 

Cells  designed  for  long  service  life,  or  high  capacity  at 
low  to  moderate  discharge  currents,  require  commensu- 
rately  high  volumes  of  active  material.  Those  designed 
for  high  discharge  rates  often  have  lower  total  volume 
of  active  material,  but  high  surface  areas.  For  example, 
the  spirally  wound,  Duracell,  CR2,  lithium  manganese 
cell,  has  a  mass  of  11  g,  volume  of  5.168  cm3,  and  ca¬ 
pacity  of  0.5  Ah  at  1  A  discharge  current,  which  has  supe¬ 
rior  higher  rate  performance  over  its  bobbin  construction 
counterpart — the  Duracell  MX1500-AA  [19]. 

(3)  environment :  temperature,  pressure,  and  exposure 
Sufficiently  low-operating  temperatures  result  in  reduced 
capacity,  and  increase  the  negative  slope  of  the  voltage 
curve  due  to  decreased  chemical  activity,  and  enhanced 
internal  resistance.  However,  high-operating  tempera¬ 
tures  can  result  in  self-discharge  [20] .  For  example,  Liaw 
et  al.  [20]  found  for  cylindrical  18650  lithium-ion  cells, 
increasing  temperature  produced  greater  capacity  and 
power  fade.  However,  as  the  state  of  charge  decreased, 
the  amount  of  power  fade  expressed  at  each  temperature 
decreased  relative  to  the  100%  state  of  charge. 
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Table  1 

Thin-film  battery  technologies  [9,10,12,13] 


Electrochemistry 

Dimensions 

Volt 

Specific  energy  (Wh  kg  1 ) 

Energy  density  (Wh  F  1 ) 

Solid  state,  thin-film  cathode  =  V2O5, 
electrolyte  =  block  polymer, 
anode  =  lithium  [13] 

Area  =  1  cm2,  thickness  =  20  pm 

1. 5-4.0 

350  @1.6C 

560  W  kg-1  @1.6C 

Thin-film,  cathode  =  LiCo02, 
anode  =  MCF 

(mesophase-pitch-based  carbon 
fiber),  electrolyte  =  organic  liquid 
[12] 

Area  =  21.7  cm2 ,  thickness  =  1.2  mm 

~3.8  V 

120  @0.2  C 

Thin-film,  cathode  =  FiCo02, 
anode  =  MCF 

(mesophase-pitch-based  carbon 
fiber),  electrolyte  =  organic  liquid 
[12] 

Area  =  21.7  cm2,  thickness  =  3.6mm 

~3.8  V 

160  @0.2  C 

Thin-film,  cathode  =  FiCo02, 
anode  =  B -doped  MCF 
(mesophase-pitch-based  carbon 
fiber),  electrolyte  =  organic  liquid 
[12] 

Area  =  21.7  cm2 ,  thickness  =1.2  mm 

~3.8  V 

172  @0.2  C 

Thin-film,  solid  state,  cathode  =  TiS2, 
anode  =  lithium, 

electrolyte  =  oxide/sulfide  glassy 
[10] 

Thickness  =  ~10  pm 

~2.4-2.5  V 

140  WhF-1  @  135  mA  cm“2 

270  WL_1 

Thin-film  rechargeable  lithium, 
cathode  =  V2O5 ,  anode  =  lithium 
metal,  electrolyte  =  amorphous 
inorganic  electrolyte  [9] 

Area  =  6.25  cm2,  cathode 
thickness  =  0.13  pm,  electrolyte 
thickness  =  1  pm,  lithium 
thickness  =  3-5  pm 

1.5-3. 5 

26.3  @.13 C 

14.4  @0.13  C 

Thin-film  rechargeable  solid  state, 
cathode  =  TiS2,  anode  =  lithium 
metal,  electrolyte  =  oxide/sulfide 
glassy  electrolyte  [10] 

Thickness  =10  pm 

1. 8-2.5  V 

140  @  75-135  |xAcm-2 

Often,  two  or  more  power  supplies  are  used  within  the 
same  system  to  overcome  some  of  the  limitations  of  elec¬ 
trochemistry,  size,  shape,  and  conditions  outlined  above. 
The  WIMS-EMT,  like  many  autonomous  MEMS  systems, 
has  low-power  standby  periods  for  sensing,  and  high-power 
pulses  for  R/F  transmission  and  reception  [24].  Typically, 
a  hybrid  power  supply  combines  a  high-power  density/low- 
energy  density  power  element  to  support  high  pulses,  and 
high-energy  density/low-power  density  element  to  provide 
low,  sustained  power  [21].  Such  systems  have  been  de¬ 
signed,  for  example,  for  automotive  applications  (e.g.  poly¬ 
mer  electrolyte  membrane  (PEMFC)  fuel  cell/lead-acid  bat¬ 
tery  studied  in  [22]),  and  civilian  power  systems  (e.g.  the 
wind/diesel/battery  system  studied  in  [23]).  These  studies, 
among  others,  have  shown  that  improved  performance  is  pos¬ 
sible  with  hybrid  supplies,  though  the  power  supplies  in  the 
cases  mentioned  were  pre-selected. 


2.  WIMS-EMT:  requirements 

The  WIMS-EMT  consists  of  a  micro-gas  chromato¬ 
graphic  (micro-GC)  system  and  MEMS  control  circuitry  de¬ 
signed  for  separation  and  analysis  of  volatile  compounds  in 
the  environment  [1].  Like  many  “wireless  devices”  presently 
under  development,  the  WIMS-EMT  uses  radio  waves  in¬ 


stead  of  current-delivering  wires,  to  transmit  signals.  The 
micro-GC  contains  two  gas  chromatographic  columns  for 
sample  separation,  one  heated  preconcentrator  for  collection 
and  concentration  of  volatile  gaseous  compounds,  a  micro¬ 
pump  and  several  valves  used  to  cycle  the  gases  through  the 
system,  a  thermal  electric  cooler  to  cool  the  sensor  array,  and 
a  sensor  array  for  the  detection  and  analysis  of  volatile  com¬ 
pounds  in  the  air.  The  MEMS  control  circuitry  contains  a 
micro-control  unit  (MCU),  a  Universal  Micro-Sensor  Inter¬ 
face  (UMSI)  chip,  and  a  radio  frequency  evaluation  board, 
all  of  which  operate  according  to  the  generic  multi-element 
MEMS  architecture,  established  as  a  platform  for  low-power 
wireless  applications  where  small  size  and  sensing  accuracy 
are  vital  [24]. 

The  power  consumption  in  the  device  spans  the  watt  3.7  W 
to  microwatt  200  jxW  ranges  (with  component  power  require¬ 
ments  listed  in  Table  2).  The  targets  for  power  system  vol¬ 
ume  and  mass  are,  respectively,  1.45  cm3  and  0.0043  kg;  the 
total  system  volume  is  ~65.08  cm3.  The  maximum  and  min¬ 
imum  (aggregate)  power  requirements  for  the  WIMS-EMT 
are  ^4.65  and  1.125  W,  respectively;  maximum  and  mini¬ 
mum  (aggregate)  discharge  currents  are  ~0.5 12  and  0.327  A. 

The  goals  of  this  paper  are  to  describe  a  systematic  method 
for  the  selection  and  design  of  power  systems  for  MEMS 
devices,  wherein  a  case  study  of  the  WIMS-ERC  will  be 
performed.  The  three  strategies  are: 
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Table  2 


WIMS-EMT  component  list  of  voltage  and  power  requirements 


Component 

Power  (mW) 

Voltage  (V) 

Chromatographic  columns  one  and  two 

450 

15 

Sensor  array 

0.20 

3 

Thermal  cooler 

200 

3 

Microvalves 

100 

1.5 

Vacuum  pump 

780 

6 

MCU 

Sustained  =  10 

0.9 

(40  MHz) 

Peak  =  22 

1.8 

Ambient  sensors 

Negligible 

3 

UMSI  chip 

12 

3 

Preconcentrator 

Sustained  =  650 

3.0 

Peak  =  3700 

16.64 

RF  evaluation  board 

Sustained  =  45 

Peak  = 125 

3 

•  specification  of  a  single,  aggregate  power  supply,  resulting 
in  a  single  battery  electrochemistry  and  cell  size; 

•  specification  of  several  power  supplies,  by  a  priori  division 
of  power  sources  by  power  range,  e.g.  the  jjlW,  mW,  and  W 
used  in  the  testbed,  resulting  in  a  hybrid  power  system  with 
a  maximum  (here)  of  three  electrochemistry  types  and  cell 
geometries;  and 

•  specification  of  an  arbitrary  number  of  power  “bundles,” 
based  on  available  space  in  the  device,  yielding  either  a 
homogenous  or  hybrid  power  system  with  one  or  more 
types  of  electrochemistry  and  cell  geometries. 

3.  Methods 

A  general  overview  of  the  algorithm  is  shown  in  Fig.  1 . 
Final  power  configurations  consist  of  electrochemical  cells 
arranged  in  series  or  parallel.  Resistors,  operational  ampli¬ 
fiers  (op-amps),  and  capacitors  are  routinely  used  to  obtain 
higher  voltage  values  or  modify  current  discharge  values  in 


Fig.  1.  Overview  of  the  power  selection  algorithm. 
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these  applications.  Other  workers  have  used  CMOS  op-amps 
[25,26],  to  alter  the  supply  voltage  for  MEMS  devices.  After 
selection  of  batteries,  we  used  specifications  of  a  set  of  three 
commercially-available  op-amps  and  one  voltage  generator 
to  reduce  the  number  of  batteries  required  only  to  drive  up 
voltage. 

A  stepwise  outline  of  the  algorithm  follows: 

1 .  User  supplies  power  system  targets  (mass,  volume,  area, 
etc.)  and  subdevice  requirements  (current  and  voltage) 
and  environmental  constraints. 

2.  Environmental  requirements  are  used  to  eliminate  elec¬ 
trochemistries. 

3.  Subdevice  performance  parameters  are  computed 
(power,  specific  energy  and  power,  and  energy  and  power 
densities). 

4.  Approaches  1,  2,  and  3  are  applied  to  subdevice  param¬ 
eters. 

5.  Energy,  voltage,  and  current  ratios  (maximum/available 
from  cell)  are  computed  for  each  battery  in  the  reduced 
database. 

6.  Conditionality  statements  are  applied  to  compute 
the  total  numbers  and  required  arrangements  of 
battery. 

7.  Power  configurations  are  eliminated  based  on  mass,  vol¬ 
ume,  or  area  target  values.  If  no  batteries  meet  these 
targets,  a  maximum  of  three  batteries  that  have  values 
closest  to  target  values  are  selected. 

8.  Batteries  are  selected  based  on  specific  or  energy  prop¬ 
erties  for  Approaches  1  and  2,  and  based  on  packaging 
constraints  for  Approach  3. 

9.  Batteries  are  selected  based  on  maximum  lifetime. 

10.  Results  from  all  three  approaches  are  compared  and  one 
system  is  selected. 

We  note  that  system  cost  was  not  used  as  a  criterion  in 
the  inherently  high-cost,  high-tech  application  studied  here, 
though  this  would  clearly  be  a  first-order  concern  in  consumer 
electronics.  The  second  approach  studied  here  was  limited 
to  evaluation  of  power  ranges  of  jjiW,  mW,  and  W  used  in 
the  testbed  for  definiteness,  though  extension  to  other  ranges 
would  be  straightforward. 


Table  3 


Input  parameters  supplied  by  the  user,  associated  variable  names,  and  data 
type 


Input  parameters 

Name 

Type 

Number  of  devices 

N 

Integer 

Current  vs.  time  for  each  device 

eft)  j'A  =N 

Real  number 

Voltage  vs.  time  for  each  device 

vj(t)j:l  =N 

Real  number 

Ambient  temperature  (K) 

T 

Real  number 

Use  of  primary  or  secondary  de- 

<p>  or  <s> 

Single  string 

vices 

character 

Volume  or  mass  optimization 

<v>  or  <m> 

Single  string 
character 

Desired  mass  for  power  supply 

T^sys 

Real  number 

(kg) 

Desired  volume  for  power  supply 

%s 

Real  number 

(cm3) 

Maximum  surface  area  for  power 

^sys 

Real  number 

supply  (cm2) 

Number  of  operation  cycles 

ncyc 

Integer 

Maximum  number  of  power  bun- 

T^loc 

Integer 

dies 

Mass  |xW  power  range 

m\ 

Real  number 

Mass  mW  power  range 

m2 

Real  number 

Mass  W  power  range 

m3 

Real  number 

Volume  |xW  power  range 

vi 

Real  number 

Volume  mW  power  range 

V2 

Real  number 

Volume  W  power  range 

V3 

Real  number 

Bundle  1  area 

a_\ 

Real  number 

Bundle  2  volume 

V2 

Real  number 

(discharge  and  charge  temperature,  pressure,  and  magnetic 
material  limitations). 

Power  requirements  per  time  increment  are  computed 
from  the  parameters  of  Table  3.  Power,  pft),  is  simply  the 
product  of  voltage  vft)  (V)  and  current,  eft)  (A),  per 

pft)  =  Ci(t)vi(t ),  i  =  1  :  A,  no  sum  (3) 

Energy,  eu  is  the  product  of  power  and  time  increment, 
given  by: 

et  =  pft)At  (4) 

specific  energy,  eu  is  energy  per  unit  mass, 

&tpiU) 

e,  = -  (5) 

mx 


3.1.  User  input  and  calculation  of  device  parameters 

A  list  of  required  input  data  is  provided  as  Table  3.  The 
algorithm  developed  here  assumed  user-provided  power  con¬ 
straints,  including  power  system  mass,  msys;  volume,  nsys; 
areal  footprint,  asys;  maximum  number  of  power  bundles, 
n\oc;  mass  range  values,  m\,  m2,  and  m3  (subscripts  1  =  jjiW, 
2  =  mW,  and  3  =  W);  and  volume  range  values,  v\ ,  V2,  and  1)3 ; 
and  mass,  volume,  and  area  bundle  values  mq,yq,  and  aq,  re¬ 
spectively,  where  q=  l:n\oc.  Also  provided  were  the  number 
of  subdevices,  A;  and  discharge  current  (cj(t))  and  voltage 
values  ( Vj(t )),  where  i  =  l:N,  for  each  time  interval  of  a  cy¬ 
cle;  battery  selection  priority  (mass  or  volume);  type  of  bat¬ 
tery  (primary  or  secondary);  and  environmental  constraints 


weighted  specific  power,  pi,  is  power  per  unit  mass, 
(A  t/tr)pi(t) 


Pi  = 


m 


(6) 


energy  density,  eu  is  energy  per  unit  volume  given  by: 
A  tiPiit) 


e;  = 


V 


(7) 


and  weighted  power  density,  pi,  is  power  per  unit  volume, 
(A  t/Pr)pi(t) 


Pi  = 


V 


(8) 


In  the  previous  expressions  tj  (s)  is  the  duration  of  a  single 
cycle.  Subscript  x-a,  corresponds  to  the  aggregate  system 
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parameters.  For  x-r ,  subscript  r  takes  values  of  1  for  the 
microwatt  (jjiW)  power  range,  2  for  the  milliwatt  (mW)  power 
range,  and  3  for  the  watt  (W)  power  range.  For  x  =  q,  we  have 
values  of  q  from  1  :n\oc ,  corresponding  to  the  number  of  power 
bundle  locations.  Eqs.  (4)— (8)  are  evaluated  independently  for 
Approaches  1,  2,  and  3,  respectively. 

3.2.  Restriction  of  electrochemistry  based  on 
environmental  constraints 

The  batteries  evaluated  by  the  algorithm  were  selected 
from  a  database  comprising  25  commercial  secondary 
(rechargeable)  batteries,  representing  the  following  elec¬ 
trochemistries:  manganese  lithium,  nickel  metal  hydride, 
lithium-ion,  and  nickel  cadmium.  One  type  of  thin-film  bat¬ 
tery  (lithium-free)  was  included  [27]  for  application  to  Ap¬ 
proaches  2  and  3.  These  commercially  available  battery 
electrochemistries  were  selected  due  to  their  relatively  low 
masses  (~4.5-73g)  and  volumes  (^0.03-56.50  cm3),  mak¬ 
ing  them  generally  suitable  for  MEMS  applications.  The  bat¬ 
tery  database  developed  for  this  work  includes  performance 
characteristics,  such  as  mass,  volume,  capacity  (at  a  defined 
discharge  rate  (C)),  nominal  discharge/charge  current,  and 
nominal  voltage,  per  Table  4. 

3.3.  Application  of  Approaches  1,  2,  and  3 

Approach  1  requires  the  summation  of  device  values  for 
power,  pi(t)  over  each  time  increment  to  produce  an  aggre¬ 
gate  device  power  versus  time  profile.  The  aggregate  power 
data  for  the  system  was  used  to  compute  aggregate  system 
parameters  such  as  energy,  Ea  and  weighted  power,  Pa,  from 
Eqs.  (9)  and  (10)  as: 

N 

Ex  =  ]T  Pj(t)t  (9) 

j=  i 


In  Approach  3,  user-defined  values  included  the  number  of 
available  power  bundle  locations,  wherein  grouping  is  based 
on  availability  of  space  and  surface  area;  the  volumes,  sur¬ 
face  areas,  and  mass  target  values  were  also  collected.  The 
required  energy  for  each  site  location  is  computed  based  on 
the  fraction  of  total  space  required  for  a  particular  bundle. 
For  example,  if  three  locations  have  volume  and  areal  foot¬ 
print  requirements  of  yi  =  2.5  cm3,  y2  -  3.5  cm3,  g3  =  5  cm3, 
=  1.0  cm2,  Q2  =  0.5  cm2,  and  a3  =  0.75  cm2,  respectively; 
then  the  fraction  of  energy  required  for  each  location  is  com¬ 
puted  as: 


V^loc  ,, 
2^q=lVq 


(ID 


where  the  energy  for  each  location  is  given  as: 


Eq  —  OiqE  a 


(12) 


where  £a  is  computed  from  Eq.  (9).  The  weighted  power 
is  computed  in  a  similar  manner,  wherein  a  fraction  of  the 
surface  areas  is  given  by: 


E"loc 

q=l-4 


(13) 


and  the  weighted  power  is  given  by: 


Pq  =  fiqPa-  (14) 

After  power  versus  time  profiles  have  been  constructed  for 
Approaches  1  and  2,  and  the  energy  and  weighted  power  have 
been  computed  for  Approach  3,  we  obtain  specific  energy  for 
Approaches  1  and  2  as: 

N 

Ex  =  ]T  e,  (15) 

i—\ 

weighted  specific  power  given  by: 


and 


N 


px(t)  =  ^2 pit* ) 
7=1 


t 

tj 


N 

Px  =  22  ^  (16) 

(10)  i=1 

energy  density  given  by: 


where  subscript  v  is  a  for  aggregate  system  of  device  param¬ 
eters;  r  for  each  power  range  (r=  1,  2,  and  3  for  jjiW,  mW, 
and  W  power  ranges,  respectively). 

In  Approach  2,  device  power  values  are  separated  accord¬ 
ing  to  power  range:  microwatt,  milliwatt,  and  watt,  wherein 
the  power  and  current  data  were  summed  over  each  time  in¬ 
crement  within  each  power  range.  This  logic  was  performed 
in  iterations:  subdevices  contributing  to  the  largest  power 
values  within  a  particular  power  range  were  removed  and 
placed  in  a  higher  power  range  than  their  initial  position  as 
needed.  Once  both  guidelines  for  the  micro  and  milliwatt 
power  ranges  are  met,  final  power  versus  time  profiles  for 
each  power  range  were  assembled  so  that  energy,  Zq,  £2,  and 
E3 ;  weighted  power,  P\,P2,  and  Pi  could  be  calculated  from 
Eqs.  (9)  and  (10). 


N 

Ex  =  yf'ei  (17> 

i=  1 

and  weighted  power  density  given  by: 

N 

*r  =  5>  US) 

i= 1 

by  summing  device  values  per  time  increment.  The  specific 
and  energy  values  for  Approach  3  are  determined  by  mul¬ 
tiplying  the  specific  energy  and  energy  density  values  com¬ 
puted  for  Approach  1  by  aq  for  each  bundle  site.  Similarly, 
the  specific  power  and  power  density  values  for  Approach  3 
can  be  calculated  by  multiplying  the  values  from  Approach 
1  by  Pq- 
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3.4.  Calculation  of  energy,  voltage,  and  current  factors 


Within  a  power  system,  total  system  energy, 


(19) 


is  based  on  placement  of  batteries;  series  placement  achieves 
higher  voltage,  and  parallel  placement  achieves  higher  dis¬ 
charge  current.  In  Eq.  (19),  bj  is  the  nominal  voltage,  with 
j=l\d,  and  where  d  is  the  number  of  batteries  selected  from 
the  database.  C/,  is  the  cell  capacity,  at  time,  t\  tj  is  the  time 
interval.  Energy  values  for  each  approach,  were  calculated  for 
the  aggregate  system,  power  ranges,  and  bundle  sites  using 
Eq.  (9). 

Cell  energy,  ej,  nominal  voltage,  bj ,  and  maximum  allow¬ 
able  discharge  current,  ij ,  are  related  to  system,  power  range, 
or  power  bundle  values  for  energy  (Ea,  Er ,  Eq),  maximum 
voltage  (Va,  Vr,  Vq),  and  maximum  discharge  current  (za,  ir, 


iq)  via: 

II 

'-*■  * 

(20) 

II 

(21) 

and 

lx 

Zj  =  T- 

li 

(22) 

All  factors  computed  from  Eqs.  (20)-(22)  are  first  rounded 
up  to  the  nearest  integer,  and  then  substituted  into  conditional 
statements. 


For  conditions  specified  in  Table  5(b),  total  number  of  bat¬ 
teries  is: 

nj  =  yjzj  +  \xj-zj\  (26) 

where  the  number  of  batteries  in  series  is  found  using  Eq.  (24), 
and  the  number  of  batteries  in  parallel  to  meet  the  discharge 
current  requirement  is  given  by: 

"'/  =  .'■/  -  './  <27) 

The  number  of  parallel-placed  batteries  required  in  to  meet 
the  energy  requirement  is  then  given  by: 

Uj  =  Xj.  (28) 

When  factor  v  is  less  than  or  equal  to  1  the  energy  provided 
from  a  single  cell  can  supply  enough  energy  for  one  cycle 
of  operation.  If  v  is  greater  than  1,  cells  must  be  placed  in 
parallel  to  meet  the  energy  requirements  while  maintaining 
operation  at  a  specific  voltage  value.  For  example,  if  v  =  3, 
then  three  batteries  placed  in  parallel  will  satisfy  the  energy 
requirement,  as  illustrated  in  Fig.  2.  Also,  if  z  is  less  than  or 
equal  to  1 ,  the  required  discharge  current  for  the  application 
can  be  met  with  one  cell,  however,  if  this  value  is  greater  than 
1,  additional  cells  placed  in  parallel  are  necessary.  Placing 
the  cells  in  parallel,  allows  the  current  draw  to  be  divided 
amongst  the  cells  to  prevent  cell  malfunction  due  to  overly 
high  discharge  rates.  Hence,  if  z  =  3,  Fig.  2  is  also  appropriate. 
If  y  is  less  than  or  equal  to  1,  then  only  1  cell  is  required  to 
satisfy  the  voltage  requirements.  However,  if  y  is  greater  than 
1,  cells  must  be  placed  in  series.  For  example,  if  y  =  3,  then 
three  cells  are  required,  as  depicted  in  Fig.  2.  In  systems 
containing  many  subdevices,  various  combinations  of  x,  y, 
and  z  can  be  used  to  satisfy  the  application. 


3.5.  Apply  conditionality  statements 


3.6.  Elimination  based  on  mass,  volume,  or  area  targets 


Conditional  statements  are  used  to  determine  the  number 
of  batteries  and  cell  configuration  required,  i.e.  placement  of 
batteries  in  series,  parallel,  and/or  a  combination  required  for 
system  operation.  Conditional  statements  have  been  derived 
for  each  possible  combination  of  x,  y,  and  z  along  with  cor¬ 
responding  circuit  diagrams  illustrated  in  Table  5(a)  and  (b) 
,  where  n a,  nr ,  and  nq  are  the  total  number  of  batteries;  sa, 
sr,  and  sq  are  the  total  number  of  batteries  in  series;  and  u;a, 
wr,  and  wq  are  the  total  number  of  batteries  in  parallel  for 
Approaches  1, 2,  and  3,  respectively.  For  conditions  specified 
in  Table  5(a),  the  total  number  of  batteries  in  the  system  is: 

nj  =  yjzj  (23) 

the  number  of  batteries  which  must  be  placed  in  series  is: 

Sj  =  yj  (24) 

and  the  number  of  batteries  which  must  be  placed  in  parallel 
is: 


(25) 


The  total  amount  of  energy  provided  by  the  system  of  cells 
given  by: 

Ej  =  einx ■  (29) 

Total  mass  and  volume,  are,  respectively, 

Mj  =  mjnx  (30) 

and 


V 


]  =  vJn 


x- 


(31) 


For  example:  for  a  system  requiring  Fa  =  0. 1  Wh,  /a  =  0.1  mA, 
na  =  3.0  V  and  for  which  two  batteries  are  examined  (Fithium, 
CR927-Renata,  with  E\  =  0.09Wh,  i\  =0.05E-3,  v\  =3.0  V; 
and  NiMH,  HHR60AAAH,  Panasonic,  with  F2  =  0.6Wh, 
?2  —  0. 1 1  A,  V2  =  1 .2  V)  we  can  calculate  the  total  number  of 
required  batteries  as  2  and  3,  respectively.  The  number  of 
cycles  a  battery  system  can  provide  is  simply  given  by: 


W1  =  Zj- 


(32) 
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Table  5 

(a)  Conditionality  statements  applicable  to  cases  where  y  or  z  has  the  maximum  value  of  the  three  factors;  (b)  conditionality  statements  applicable  to  cases 
where  x  has  the  maximum  value  of  the  three  factors,  or  where  x  is  greater  than  z 


Condition  Expression 

(a) 


x  =  y  =  z 


x<y  <z  nj=yz,  Sj=y,  wj=z 


y<x<z 


x<z<y 


x  =  y<z 


x  =  z<y 


x  =  z>y 


y  =  z>x 


Examples  and  circuit  diagram  x  y  z  nj  sj  Wj 
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Table  5  ( Continued ) 


Condition  Expression 

(b) 

z  <  x  <  y 


rij  =yz+\x  -  z |;  sj=y; 

Wj=y  —  z;  iij  =  x 
x=y>z  ' 


y  <  z  <  x 


z<y  <x 


y  =  z<x 


Examples  and  circuit  diagram 


x 


y 


nl 


h 

& 


Kt>C\2> 


J-EH 


3 

2 


5 

2 


2 

1 


11 

3 

7 


7 


P 


|cf> 


Si 

3 

2 

2 

4 

1 

2 

2 


Wj  and  Uj 

1  and  2 

2  and  3 

1  and  2 

2  and  3 

0  and  3 

3  and  4 

1  and  3 


3  2  and  5 

1  0  and  3 

2  2  and  4 


Battery  configurations  that  do  not  meet  the  mass,  msys,  vol¬ 
ume,  nsys,  or  areal  footprint,  asys  targets  are  eliminated,  by 
our  algorithm,  from  the  list  of  potential  cell  configurations.  If 
no  batteries  remain  after  this  elimination  step,  battery  config¬ 
urations  having  the  least  mass,  volume,  or  areal  (depending 
on  user  specification)  values  are  selected. 


3. 7.  Selection  based  on  specific  energy  or  energy  density 

If  volume  prioritization  is  selected  by  the  user,  energy 
density  is  computed  using  the  expression: 

Ej  =  *5-  (33) 

Vj 


Fig.  2.  Circuit  diagrams  illustrating  batteries  in  both  parallel  and  series. 
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for  the  remaining  batteries.  Cell  configurations  providing  the 
maximum  energy  density  are  selected  next.  If  mass  prioritiza¬ 
tion  is  selected,  battery  selection  is  based  on  specific  energy 
as: 


where  the  maximum  specific  energy  generated  from  this  ex¬ 
pression  corresponded  to  battery(s)  selected. 

3.8.  Selection  based  on  battery  lifetime 

If  more  than  one  battery  remains  after  application  of  all 
previous  constraints,  then  they  were  subjected  to  evaluation 
of  lifetime.  Configurations  providing  the  longest  battery  life¬ 
time  for  a  particular  approach  are  used.  This  lifetime  is: 

Lj  =  tMM(Nj)  (35) 

where  p  equal  to  1  for  primary  cells  or  equal  to  the  number 
of  cycles  a  secondary  battery  can  provide.  The  final  results 
for  each  approach  are  compared  and  an  optimal  solution  is 
selected. 

3.9.  WIMS/ERC  testbed 

Power  and  voltage  requirements  for  components  of  the 
WIMS-EMT  are  listed  in  Table  2.  User  input  parameters  are 
listed  in  Table  6.  Since  no  environmental  constraints  were 
provided  for  the  WIMS-EMT  for  discharge  and  charge  tem¬ 
perature,  charge  rate,  method  of  charge,  storage  time,  etc., 
all  battery  electrochemistries  were  considered  potential  can¬ 
didates,  and  the  algorithm  analysis  began  with  calculation  of 
device  parameters. 

4.  Results 

4.1.  Approach  1 

4.1.1.  Calculated  device  parameters 

Energy  and  power  data  for  each  device  was  summed  per 
time  segment  to  produce  aggregate  system  parameters  for 
Approach  1.  These  values  for  energy  (0.39  Wh),  weighted 


Table  6 


EMT  input  constants  and  power  supply  design  targets 


Number  of  device  in  the  systems,  N 

10 

Ambient  temperature  [°C] 

20 

[p] rimary  or  [secondary  battery  system 

S 

[vjolume  or  [m]ass  optimization 

M 

Desired  mass  [kg] 

0.0043 

Desired  volume  [cm3] 

1.45 

Number  of  cycles 

1 

Number  of  power  bundles,  n\oc 

4 

Mass  [kg]  p,W  range,  m\ 

0.000086 

Mass  [kg]  mW  range,  m2 

0.00077 

Mass  W  range,  m3 

0.00344 

Volume  [cm3]  |xW  range,  v\ 

0.029 

Volume  [cm3]  mW  range,  V2 

0.261 

Volume  [cm3]  W  range,  V3 

0.00344 

Volume  [cm3]-bundle  1 

3.68 

Area  [cm2] — bundle  2 

61.0 

power  (1.42  W),  specific  energy  (91.06  Wh  kg-1),  energy 
density  (270.04  WhL-1),  specific  power  (330.49  W  kg-1), 
and  power  density  (980.0  W  L- 1 )  system  are  listed  in  Table  7. 

4.1.2.  Calculation  of  energy,  voltage,  and  current 
factors — application  of  conditionality  statements 

For  each  of  the  25  batteries  in  the  database,  energy  {xf), 
voltage  (yj),  and  current  (zj)  factors  were  computed;  along 
with  the  total  number  of  batteries  required,  nj ,  number  in 
series,  Sj,  number  in  parallel,  Wj ,  configuration  mass,  Mj, 

/V  — 

volume,  Vj,  number  of  cycles,  Nj,  and  lifetime,  Lj.  The  results 
from  these  calculations  are  listed  in  Table  8. 

4.1.3.  Selection  based  on  mass  requirements 

The  WIMS-EMT  requires  mass  prioritization,  where 
the  target  mass  of  the  power  system  is  0.0043  kg.  Clearly, 
from  the  results  listed  in  Table  8,  no  battery  system  in  the 
database  met  this  target  value,  so  the  systems  producing 
the  minimum  mass  values  of  0.068  and  0.069  kg  (Ultral¬ 
ife  UBC64 1730/PCM  and  UBP3 83450/PCM)  were  selected. 
For  the  Ultralife  UBC64 1730/PCM,  15  batteries  were  re¬ 
quired  ( n  =  15),  with  5  in  series  (s  =  5),  in  3  parallel  rows 
( w  =  3).  Selection  of  the  Ultralife  UBP383450/PCM  config¬ 
uration  required  n  =  5;  all  were  placed  in  series  (s  =  5).  A 
basic  circuit  diagram  of  the  power  system  design  based  on 
Approach  1  is  detailed  in  Fig.  3. 


Table  7 


Approach  1 :  power  and  energy  values  as  a  function  of  time 


Time 

Power 

(Watts) 

Calculated 
voltage  (V) 

Current 

(A) 

Energy 

(Wh) 

Specific 
energy 
(Wh  kg- 1 ) 

Energy 

density 

(WhL-1) 

Weighted 
power  (W) 

Specific 

power 

(Wkg-1) 

Power 

density 

(WL->) 

Maximum 
voltage  (V) 

5 

1.125 

3.438 

0.327 

0.002 

0.363 

1.077 

0.006 

1.319 

3.911 

6 

555 

1.115 

3.527 

0.316 

0.172 

39.959 

118.499 

0.624 

145.073 

430.218 

6 

3 

4.735 

9.252 

0.512 

0.004 

0.918 

2.722 

0.014 

3.330 

9.876 

17 

9 

1.91 

5.56 

0.344 

0.005 

1.112 

3.298 

0.017 

4.030 

11.951 

15 

420 

1.795 

3.869 

0.464 

0.209 

48.707 

144.443 

0.760 

176.739 

524.124 

15 

992 

10.68 

Total 

0.392 

91.060 

270.039 

1.421 

330.492 

980.079 

Mass  (kg)  =  0.0043;  volume  (L)  =  0.00145. 
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Fig.  3.  Circuit  diagram  for  solution  proposed  from  Approach  1. 


4.1.4.  Selection  based  on  specific  energy 

Of  the  two  battery  systems  examined,  the  Ultral¬ 
ife  UBC64 1730/PCM  had  the  highest  specific  energy, 
164.44  Wh  kg-1,  although  its  predicted  lifetime  (2.4E03h) 
was  much  less  than  a  system  designed  using  the 
UBP3 83450/PCM  cell  (4.1E03  h). 

4.1.5.  Selection  based  on  battery  lifetime 

Ultralife  UBP3 8 3450/PCM  battery  provided  a  longer  life 
(4.1E03  h)  in  comparison  to  UBC641730/PCM  (2.4E03h). 
The  resulting  total  number  of  batteries  needed  is  15,  and  the 
total  mass  and  volume  of  the  power  system  is  0.068  kg  and 
0.058  L,  respectively. 

4.2.  Approach  2 

4.2.1.  Calculated  system  parameters 

The  WIMS-EMT  device  power  and  energy  values  are 
summarized  in  Table  9.  Arrangement  #2  was  chosen  so 
that  the  power  range  for  the  mW  section  was  less  than 
1  W;  subdevices  with  the  largest  power  values  in  the  mil¬ 
liwatt  range  were  moved  from  the  milliwatt  power  range  to 
the  watt  power  range.  Energy,  power,  specific  energy,  and 
energy  density  values  for  each  power  range  are  shown  in 
Table  10.  Energy,  power,  specific  energy,  and  energy  den¬ 
sity  values  for  each  power  range  are  shown  in  Table  10. 
Resulting  energy  and  power  values  for  the  |xW  range  were 
0.023  ptWh  and  0.085  mW,  respectively,  and  the  specific  en¬ 
ergy  and  energy  density  values  were  0.271  Wh  kg-1  and 
0.803  Wh  L_1 ,  respectively.  The  energy  and  power  values  for 
the  mW  range  were  0.13  Wh  and  0.47  W,  respectively,  and  the 
specific  energy  and  energy  density  values  168.048  Wh  kg-1 
and  498.350  Wh  L~ 1 .  Lastly,  energy  and  weighted  power  val¬ 


ues  for  the  watt  power  range  were  0.27  Wh  and  0.99  W,  re¬ 
spectively,  and  the  specific  energy  and  energy  density  values 
were  79.275  Wh  kg"1  and  287.690  WhL-1. 

4.2.2.  Calculation  of  energy,  voltage,  and  current 
factors — application  of  conditionality  statements 

Research  on  development  of  high-power  and  capacity 
thin  film  and  microbattery  technologies  has  intensified  in 
recent  years  (e.g.  [9,10,27,18]).  We  studied  implementation 
of  lithium- free  thin  film  batteries  proposed  by  Neudecker 
et  al.  [27],  for  the  power  source  for  the  micro- watt  power 
range.  According  to  Neudecker  et  al.  [27],  these  batteries 
have  the  capability  of  cycling  between  4.2  and  3.0  V,  and  have 
demonstrated  1000  cycles  at  ImA/cm2,  or  over  500  cycles 
at  5  mA/cm2.  Also  when  operated  below  0.1  mA/cm2,  they 
provided  a  capacity  of  about  200  mAh/cm2  from  a  Li-free 
battery  with  3.2  mm  thick  LiCo02  cathode.  Hence,  this  bat¬ 
tery  technology  was  added  to  the  25  commercial  batteries  in 
the  database  for  the  power  selection  for  the  jxW  power  range. 
The  energy  (xj),  voltage  (yj),  and  current  ( Zj )  factors  were 
computed;  along  with  the  total  number  of  batteries  required, 
nj  number  in  series,  sj,  number  in  parallel,  Wj,  configuration 

A  — 

mass,  Mj,  volume,  V),  number  of  cycles,  Nj,  and  lifetime,  Lj 
for  the  all  power  ranges.  The  results  for  the  jxW,  mW,  and  W 
ranges  are  detailed  in  Table  ll(a-c),  respectively. 

4.2.3.  Selection  based  on  mass  requirements 

The  thin-film  technology  is  the  obvious  choice  for 
the  jxW  power  range,  as  its  mass  and  volume  are  very 
small,  and  its  surface  area  is  1  cm2.  A  lithium  polymer 
cells,  UBC422030/PCM  and  UBC641730/PCM  provided  the 
smallest  power  system  mass  values:  0.008  and  0.009  kg  for 
the  mW  power  range.  Also  two  lithium  polymer  batteries, 
UBC64 1730/PCM  and  UBC422030/PCM  provided  the  low¬ 
est  mass  results  for  the  W  power  range  as  well,  0.023  and 
0.040  kg,  respectively. 

4.2.4.  Selection  based  on  specific  energy  requirements 

Of  the  two  batteries  proposed  for  the  mW  power  range, 
the  UBC64 1730/PCM  battery  is  selected  because  its  specific 
energy,  164.44  Wh  kg-1  is  greater  than  the  specific  energy  of 


Table  9 

Arrangements  #1  and  #2  for  Approach  2 

Arrangement  #1  Arrangement  #2 


Power 

range 

Device  list 

Power 

(mW) 

Total  energy 
(mWh) 

Power 

range 

Device  list 

Power 

(mW) 

Total  energy 
(mWh) 

|xW 

GC  sensor  array 

0.085 

0.023 

|xW 

GC  sensor  array 

0.085 

0.023 

mW 

TE  cooler,  micro  valves,  MCU, 
vacuum  pump,  RF  board,  col¬ 
umn  #1 ,  column  #2,  preconcen¬ 
trator 

1294 

357 

mW 

TE  cooler,  microvalves  MCU, 
RF  board 

472 

130 

W 

Preconcentrator 

11 

3 

W 

Column  #1 ,  column  #2,  precon¬ 
centrator,  vacuum  pump 

990 

273 

194 
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Table  10 


Approach  2:  micro,  milli,  and  Watt  range  power  and  energy  values  as  a  function  of  time 


Time 

Microwatt 

Microwatt 

Calculated 

Energy  (Wh) 

Weighted 

Specific 

Energy 

Specific 

Energy 

range 

range  (W) 

voltage  (V) 

power  (W) 

energy 

density 

power 

density 

current  (A) 

(Wh  kg- 1 ) 

(Wh  L- 1 ) 

(Wkg-1) 

(WhL-1) 

Microwatt  range 

613 

0 

0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

419 

6.67E-05 

2.00E-04 

3.0 

2.33E-05 

0.0001 

0.27 

0.80 

0.98 

2.91 

Time 

Discharge 

Milliwatt 

Energy 

Weighted 

Specific 

Energy 

Specific 

Power 

current  (A) 

power  (W) 

(Wh) 

power  (W) 

energy 

density 

power 

density 

(Wh  kg-1) 

(Wh  L- 1 ) 

(Wkg-1) 

(WL->) 

Milliwatt 

rangea 

5 

0.0022 

0.3470 

0.0005 

0.0017 

168.05 

498.35 

609.85 

1808.53 

555 

0.12795 

0.34500 

0.0532 

0.1930 

3 

0.12795 

0.34500 

0.0003 

0.0010 

9 

0.12795 

0.34500 

0.0009 

0.0031 

420 

0.23469 

0.64500 

0.0753 

0.2731 

Total 

0.13 

0.47 

Time 

Discharge 

Milliwatt 

Energy 

Weighted 

Specific 

Energy 

Specific 

Power 

current  (A) 

power  (W) 

(Wh) 

power  (W) 

energy 

density 

power 

density 

(Wh  kg-1) 

(Wh  L“ 1 ) 

(Wkg-1) 

(WL->) 

Watt  rangeb 

5 

0.12987 

0.78000 

0.0011 

0.0039 

79.27 

235.09 

287.69 

853.15 

555 

0.12987 

0.78000 

0.1203 

0.4364 

3 

0.13051 

4.48000 

0.0037 

0.0135 

9 

0.13652 

1.65500 

0.0041 

0.0150 

420 

0.14216 

1.23000 

0.1435 

0.5208 

Total 

0.27 

0.99 

a  m2  (kg)  =  0.000774,  v2  (L)  =  0.0003,  V2  (V)  =  3.0. 
b  mx  (kg)  =  0.00344,  vi  (L)  =  0.0012,  Vi  (V)=17.0. 


UBC422030/PCM,  138.75.  Also,  battery  UBC641730/PCM 
provided  sufficiently  greater  lifetime  (979  h)  compared  to 
UBC422030/PCM  (734  h).  For  mW  power  range,  where 
UBC641730/PCM  is  used,  ^2  =  2,  S2  =  l,  and  u>2  =  2,  indi¬ 
cating  that  a  total  of  two  batteries  are  required  to  be  placed 
in  parallel  for  the  mW  power  range.  Of  the  two  batteries 
proposed  for  the  W  power  range,  UBC422030/PCM  and 
UBC64 1730/PCM,  the  latter  provides  the  largest  specific  en¬ 
ergy,  164Whkg_1  compared  to  139Whkg_1  provided  by 
UBC422030/PCM.  The  UBC422030/PCM,  however,  does 
provide  a  longer  battery  life  (1.75E03h)  in  comparison  to 
UBC64 1730/PCM  (1. 17E03  h).  For  this  power  range,  ^3  =  5, 
S3  =  5,  and  m3  =  1 ,  wherein  a  total  of  five  batteries  are  required 
to  be  placed  in  series. 

4.2.5.  Selection  based  on  battery  lifetime 

All  battery  electrochemistry  and  type  selection  were  deter¬ 
mined  in  Section  4.2.4,  however,  since  this  is  a  hybrid  system, 
battery  lifetimes  vary  depending  on  the  power  range.  The  bat¬ 
tery  lifetime  predicted  for  the  jjiW,  mW,  and  W  power  ranges 
are  8.89E03,  979,  and  164  h.  The  total  number  of  batteries 
required  is  eight,  wherein  the  total  mass  and  volume  of  the 
power  system  based  on  Approach  2  is  0.032  kg  and  0.028  L, 
respectively. 


4.3.  Approach  3 

4.3.1.  Calculated  parameters 

The  largest  surface  area  available  for  application  of  bat¬ 
teries  on  the  WIMS-EMT  is  60.1  cm2,  along  the  bottom  of 
the  PC  Board.  Since  packaging  is  still  under  development, 
we  use  this  area  to  define  one  power  bundle,  and  desig¬ 
nate  one  other  location  for  a  high-energy  density  battery 
that  can  potentially  be  stacked  on  top  of  the  device  or  along 
side  of  the  device,  where  the  volume  will  be  limited  to  1/4 
(3.7  cm3)  the  size  of  the  fluidic  substrate  which  spans  the 
entire  length  of  the  device.  The  strategy  is  to  use  a  thin- 
film  battery  to  assist  in  providing  current  for  the  high-power 
pulses  of  the  power  versus  time  curve  and  to  use  a  high- 
energy  density  cell  to  provide  the  moderate  current  loads. 
A  schematic  demonstrating  the  concept  of  grouping  power 
supplies  by  location  is  depicted  in  Fig.  4.  The  amount  of  en¬ 
ergy  required  to  supply  bundle  1  is  0.38  Wh  and  the  weighted 
power  is  1.38  W  (Table  12).  Hence,  estimates  for  the  energy 
and  power  densities  are  103.6  and  376.6  WL-1  for  bundle 
1 ,  respectively.  The  energy  and  weighted  power  for  bundle  1 
are  0.08  and  0.17  W,  respectively.  The  energy  per  unit  area 
and  power  per  unit  area  required  for  bundle  2  are  8E-04  and 
2.8E-03  Wcm-2. 


Table  1 1 

Results  from  application  of  Approach  2  (a)  microwatt  (b)  milliwatt  (c)  watt  power  ranges 


Microwatt  Power  Range 


Electrochemistry 

Part  Number 

Discharge 
Current  (A) 

Energy 

[Wh] 

Nominal 
Voltage  (V) 

X 

y 

Z 

n 

S 

w 

M  (kg) 

V  (L) 

area(cm2) 

Number  of 
Cycles  (not 
including 
recharge 
cycles) 

Lifetime 

Energy 

Density 

[Wh/cnf] 

Lithium  Fee  Thin-Film  Battery 

- 

1.00E-04 

7  2E-01 

3.60 

0.00 

0.83 

0.67 

i 

1 

i 

- 

- 

i 

30930  79 

8.87E+06 

0720 

Electrochemistry 

Part  Number 

Discharge 
Current  (A) 

Energy 

[Wh] 

Nominal 
Voltage  (V) 

X  (energy 
ratio) 

y 

(voltage 

ratio) 

Z  (current 
ratio) 

n  (number 
of  batteries) 

S  (series) 

w 

(parallel) 

M  (kg) 

Number  of  Cycles 
(not  including  recharge 
cycles) 

Lifetime  [h] 

mu 

Manganese  Lithium 

ML414S 

5.00E-06 

0.004 

3.0 

6.47E-03 

1.0 

13.3 

14 

1.0 

14  0 

0.001 

3.55E-04 

2165.2 

6.21  E+05 

45.00 

142  10 

ML421S 

3.00E-06 

0.007 

3.0 

3.37E-03 

1.0 

222 

23 

1.0 

23.0 

0.003 

8.74E-04 

6817.7 

1  95E+06 

62  73 

181  58 

HHR60AAAH 

1.00E-01 

0.600 

1.2 

3.88E-05 

2.5 

6.67E-04 

3 

3.0 

1.0 

0.036 

1.16E-02 

77327.0 

1.11E+07 

50.00 

155  71 

HHR70AAAJ 

1.44E-01 

0.864 

1.2 

2.69E-05 

2.5 

4.63E-04 

3 

3.0 

1.0 

0.039 

1.16E-02 

111350  8 

1 60E+07 

66.46 

224  23 

NH50-D 

500E-01 

3  000 

1.2 

776E-06 

25 

1  33E-04 

3 

3.0 

1.0 

0219 

1  69E-01 

386534  8 

554E+07 

41  10 

53  10 

NH35-C 

5.00E-01 

3.000 

1.2 

7.76E-06 

2.5 

1.33E-04 

3 

3.0 

1.0 

0.180 

8  09E-02 

386634  8 

5.54E+07 

50.00 

111  29 

Nickel  Metal  Hydride 

NH12-AAA 

1.70E-01 

1.020 

1.2 

2.28E-05 

2.5 

3.92E-04 

3 

mm 

1.0 

0.036 

1.16E-02 

131455  8 

1.88E+07 

85.00 

264  71 

NH15-AA 

460E-01 

2  760 

1.2 

8  43E-06 

2.5 

1  45E-04 

3 

3.0 

1.0 

0081 

2  50E-02 

355704  1 

5.10E+07 

102  22 

330  97 

NH22-9V 

3.00E-02 

0.180 

1.2 

1  29E-04 

2.5 

2  22E  03 

3 

mm 

1.0 

0.123 

6.52E-02 

23198  1 

3.33E+06 

4.39 

8.29 

SP75AAAMH 

1.50E-01 

0.900 

1.2 

2.59E-05 

2.5 

444E-04 

3 

mm 

1.0 

0.039 

1  OOE-02 

115990  5 

1  66E+07 

69  23 

269  63 

SP180AAKH 

3.60E-01 

2.160 

1.2 

1.08E-05 

2.5 

1.85E-04 

3 

mm 

1.0 

0.084 

2  38E-02 

278377  1 

3.99E+07 

77.14 

272  78 

SP230AAKH 

4.60E-01 

2  760 

1.2 

8.43E-06 

2.5 

1.45E-04 

3 

3.0 

1.0 

0.090 

2  31E-02 

355704  1 

1.02E+08 

92.00 

358  59 

CGR 17500 

1.56E-01 

2  988 

3.6 

7.79E-06 

0.8 

4.27E-04 

1 

1.0 

1.0 

0.025 

1.11E-02 

128362  8 

1 84E+07 

119  52 

268  56 

CGR18650HG 

3.40E-01 

6.480 

3.6 

3  59E-06 

0.8 

1.96E-04 

1 

1.0 

1.0 

0.042 

1.75E-02 

278377  1 

3  99E+07 

154  29 

370  88 

UBC64 1730/PCM 

1.50E-01 

0.740 

3.7 

3.15E-05 

0.8 

4  44E-04 

1 

1.0 

1.0 

0.005 

3  85  E- 03 

31790.0 

2.73E+06 

164  44 

19220 

UBP3834  50/PCM 

300E-01 

2  220 

3.7 

1  05E-05 

08 

222E-04 

1 

1.0 

10 

0  014 

8  55E-03 

95369  9 

1  37E+07 

160  87 

259  54 

Lithium  Ion 

UBC422030/PCM 

7.50E-02 

0.555 

3.7 

4.19E-05 

08 

8  89E-04 

1 

1.0 

1.0 

0.004 

3.06E-03 

23842  5 

2.05E+06 

13875 

181  39 

ICP383450G 

1.40E-01 

2  590 

3.7 

8.99E-06 

0.8 

4.76E-04 

1 

1.0 

1.0 

0.015 

6.59E-03 

111264  9 

1.59E+07 

17862 

39301 

ICP043048G 

1  44E-01 

2.564 

3.7 

8.74E-06 

0.8 

463E-04 

1 

1.0 

1.0 

0.015 

6.49E-03 

114443  9 

1 64E+07 

183  72 

410  70 

CGA523450A 

0.9 

3.3 

3.6 

7.1  IE-06 

0.8 

7.41  E-05 

1 

1.0 

1.0 

0.020 

8.93E-03 

140735  1 

2.02E+07 

168  00 

367  06 

CGA523436 

0.7 

2.5 

3.6 

9.24E-06 

0.8 

9.80E-05 

1 

1.0 

1.0 

0.015 

6.36E-03 

108257  8 

1  55E+07 

173  79 

395  93 

P-11AAH 

2.40E-02 

0.144 

1.2 

1.62E-04 

2.5 

2.78E-03 

3 

KB 

1.0 

0.020 

8.67E-03 

18558.5 

2.66E+06 

22.15 

49.83 

Nickel  Cadmium 

P-18N 

3.80E-02 

0  228 

1.2 

1.02E-04 

2.5 

1.75E-03 

3 

mm 

1.0 

0.024 

1  02E-02 

29384  2 

4.21E+06 

28.50 

67.20 

SP30AA2/3MC 

1  50E+00 

0  360 

1.2 

6  47E-05 

2.5 

4  44E-05 

3 

30 

1.0 

0033 

1  35E-02 

46396  2 

6  65E+06 

3273 

7976 

SP28AAAMC 

1  40E+00 

0.336 

1.2 

6.S3E-05 

2.5 

4.76E-05 

3 

mm 

1.0 

0.027 

1  04E-02 

43303  1 

6.21  E+06 

37.33 

96.98 
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Fig.  4.  Schematic  illustrating  implementation  of  Approach  3,  grouping  by  power  site  location. 


Table  12 


Approach  3:  power  and  energy  data  for  power  bundle  sites 


Time 

Power 

Calculated 

Current 

Energy 

Weighted 

Maximum 

Power 

Energy 

(Watts) 

voltage 

(A) 

(Wh) 

power  (W) 

voltage  (V) 

density 

density 

(V) 

(W  L-1) 

(WhL-1) 

Power  bundle  #la 

5 

1.089846 

3.438 

0.317 

0.002 

0.005 

6 

376.571 

103.766 

555 

1.079262 

3.527 

0.306 

0.166 

0.604 

6 

3 

4.644504 

9.252 

0.502 

0.004 

0.014 

17 

9 

1.85704 

5.56 

0.334 

0.005 

0.017 

15 

420 

1.756526 

3.869 

0.454 

0.205 

0.744 

15 

992 

Total 

0.381 

1.384 

Bundle  2  current  (A) 

Bundle  2 

Energy 

Weighted 

Energy 

Power 

power  (W) 

(Wh) 

power  (W) 

density 

density 

(Wh  cm-2) 

(W  cm-2) 

Power  bundle  #2b 

0.06 

1.02 

2.92E-01 

1.02 

4.79E  -  03 

1.67E-2 

a  vi  (L)  =  3.68E  -  03,  V\  (V)  =  17. 
b  a2  (cm2)  =  6.10E  +  01,  V2  (V)=17.0. 


4.3.2.  Calculation  of  energy,  voltage,  and  current 
factor s-application  of  conditionality  statements 

The  maximum  number  of  thin  films  available  for  use 
is  60  wherein  each  film  is  approximately  1  cm2  and  is 
capable  of  providing  5mA/cm2.  The  energy  (xj),  volt¬ 
age  (yj),  and  current  (zj)  factors  were  computed;  along 
with  the  total  number  of  batteries  required,  nj ,  num¬ 
ber  in  series,  sp  number  in  parallel,  Wj ,  configuration 

/V  — 

mass,  Mj ,  volume,  Vj,  number  of  cycles,  Nj ,  and  life¬ 
time,  Lj  for  all  the  power  bundle  locations  and  the  re¬ 
sults  are  in  Table  13(a  and  b)  for  bundle  1  and  2, 
respectively. 

4.3.3.  Selection  based  on  mass  requirements 

Results  for  Approach  3  are  listed  in  Table  13(a) 
bundle  1  and  (b)  bundle  2.  The  minimum  masses 
for  bundle  1  are  obtained  from  UBC64 1730/PCM 
(0.068  kg),  UBP383450/PCM  (0.069  kg),  and  CGA523436 
(0.073  kg).  Thin-film  lithium-free  batteries  are  used  for 
bundle  2. 


4.3.4.  Selection  based  on  specific  energy  requirements 
The  Panasonic  battery,  CGA523436  (0.073  kg),  was  se¬ 
lected  for  two  reasons.  First,  it  had  the  highest  specific  energy 
(173.70  Wh  kg-1)  of  the  all  three  batteries,  and  the  smallest 
volume. 

4.3.5.  Selection  based  on  battery  lifetime 

All  battery  electrochemistry  and  type  selection  are  com¬ 
pleted  in  Section  4.3.3,  however,  since  this  is  a  hybrid  system, 
battery  lifetimes  vary  depending  on  the  power  bundle.  The 
battery  lifetime  predicted  for  bundles  1  and  2  are  2. 1E04  and 
4.74E03  h,  respectively.  The  resulting  total  number  of  batter¬ 
ies  required  is  65,  and  the  total  mass  and  volume  of  the  power 
system  is  0.073  kg  and  0.032  L. 

5.  Implementation  of  operational  amplifier  (op-amp) 
and  voltage  generator  technology 

The  methodology  detailed  here  assumes  no  implementa¬ 
tion  of  additional  devices  (e.g.  resistors,  voltage  generators  or 


Table  13 

Power  bundle  results  for  (a)  bundle  1  and  (b)  bundle  2 


Electrochemi  stry 

Part  number 

Discharge 
current  (A) 

Energy 

(Wh) 

Nominal 
voltage  (Y) 

x  (energy 
ratio) 

y  (voltage 
ratio) 

z  (current 
ratio) 

n  (number  of 
batteries) 

5  (series) 

w  (parallel)  M  (kg)  Vj  (L) 

Number  of  cycles  (not 
including  recharge 
cycles) 

Lifetime  (h) 

Specific 

energy 

(Whkg-1) 

Energy  density 
(WhL-1) 

Bundle  1 

Man¬ 

ML414S 

5.00E  - 

06 

0.004 

3.0 

105.9 

5.7 

88400.00 

530400.00 

6.00 

88400.00 

42.432 

13.4370 

5.01E  +  03 

1.44E  +  06 

45.00 

142.10 

ganese 

lithium 

ML421S 

3.00E  - 

06 

0.007 

3.0 

55.3 

5.7 

147333.33 

884004.00 

6.00 

147334.00 

97.240 

33.593 

1.60E  +  04 

4.59E  +  06 

62.73 

181.58 

HHR60AAAH 

1.00E  - 

01 

0.600 

1.2 

0.6 

14.2 

4.42 

75.00 

15.00 

5.00 

0.900 

0.289 

118.005 

1.69E  +  04 

50.00 

155.71 

HHR70AAAJ 

1.44E  - 

01 

0.864 

1.2 

0.4 

14.2 

3.07 

60.00 

15.00 

4.00 

0.780 

0.231 

135.941 

1.95E  +  04 

66.46 

224.23 

NH50-D 

5.00E  - 

01 

3.000 

1.2 

0.1 

14.2 

0.88 

15.00 

15.00 

1.00 

1.095 

0.847 

118.005 

1.69E  +  04 

41.10 

53.10 

Nickel 

NH35-C 

5.00E  - 

01 

3.000 

1.2 

0.1 

14.2 

0.88 

15.00 

15.00 

1.00 

0.900 

0.404 

118.005 

1.69E  +  04 

50.00 

111.29 

metal 

NH12-AAA 

1.70E  - 

01 

1.020 

1.2 

0.4 

14.2 

2.60 

45.00 

15.00 

3.00 

0.540 

0.173 

120.365 

1.73E  +  04 

85.00 

264.71 

hy¬ 

NH15-AA 

4.60E  - 

01 

2.760 

1.2 

0.1 

14.2 

0.96 

15.00 

15.00 

1.00 

0.405 

0.125 

108.564 

1.56E  +  04 

102.22 

330.97 

dride 

NH22-9V 

3.00E  - 

01 

0.180 

1.2 

2.1 

14.2 

14.73 

225.00 

15.00 

15.00 

9.225 

4.887 

106.204 

1.52E  +  04 

4.39 

8.29 

SP180AAKH 

1.50E  - 

01 

0.900 

1.2 

0.4 

14.2 

2.95 

45.00 

15.00 

3.00 

0.585 

0.150 

1.06E  +  02 

1.52E  +  04 

69.23 

269.63 

SP230AAKH 

3.60E  - 

01 

2.160 

1.2 

0.1 

14.2 

1.23 

30.00 

15.00 

2.00 

0.840 

0.238 

1.70E  +  02 

2.44E  +  04 

77.14 

272.78 

SP230AAKH 

4.60E  - 

01 

2.760 

1.2 

0.1 

14.2 

0.96 

15.00 

15.00 

1.00 

0.450 

0.115 

1.09E  +  02 

3.11E  +  04 

92.00 

358.59 

CGR17500 

1.56E  - 

01 

2.988 

3.6 

0.1 

4.7 

2.83 

15.00 

5.00 

3.00 

0.375 

0.167 

1.18E  +  02 

1.68E  +  04 

119.52 

268.56 

CGR18650HG 

3.40E  - 

01 

6.660 

3.7 

0.1 

4.6 

1.30 

10.00 

5.00 

2.00 

0.420 

0.175 

1.75E  +  02 

2.50E  +  04 

158.57 

381.18 

UBC64 1730/PCM 

1.90E  - 

01 

0.740 

3.7 

0.5 

4.6 

2.33 

15.00 

5.00 

3.00 

0.068 

0.058 

2.91E  +  01 

2.50E  +  03 

164.44 

192.20 

Lithium 

UBP383450/PCM 

5.90E  - 

01 

2.220 

3.7 

0.2 

4.6 

0.75 

5.00 

5.00 

1.00 

0.069 

0.043 

2.91E  +  01 

4.17E  +  03 

160.87 

259.54 

ion 

UB  C422030/PCM 

1.40E  - 

01 

0.555 

3.7 

0.7 

4.6 

3.16 

20.00 

5.00 

4.00 

0.080 

0.061 

2.91E  +  01 

2.50E  +  04 

138.75 

181.39 

ICP383450G 

1.40E  - 

01 

2.590 

3.7 

0.1 

4.6 

3.16 

20.00 

5.00 

4.00 

0.290 

0.132 

1.36E  +  02 

1.95E  +  04 

178.62 

393.01 

ICP383450G 

1.44E  - 

01 

2.664 

3.7 

0.1 

4.6 

3.07 

20.00 

5.00 

4.00 

0.290 

0.130 

1.40E  +  02 

2.00E  +  04 

183.72 

410.70 

CGA523450A 

9.00E  - 

01 

3.276 

3.6 

0.1 

4.7 

0.49 

5.00 

5.00 

1.00 

0.098 

0.045 

4.30E  +  01 

6.16E  +  03 

168.00 

367.06 

CGA523436 

6.80E  - 

01 

2.520 

3.6 

0.2 

4.7 

0.65 

5.00 

5.00 

1.00 

0.073 

0.032 

3.30E  +  01 

4.74E  +  03 

173.79 

395.93 

Nickel 

P-11AAH 

2.40E  - 

02 

0.144 

1.2 

2.6 

14.2 

18.42 

285.00 

15.00 

19.00 

1.853 

0.824 

1.08E  +  02 

1.54E  +  04 

22.15 

49.83 

cad¬ 

P-18N 

3.80E  - 

02 

0.228 

1.2 

1.7 

14.2 

11.63 

180.00 

15.00 

12.00 

1.440 

0.611 

1.08E  +  02 

1.54E  +  04 

28.50 

67.20 

mium 

SP30AA2/3MC 

1.50E  +  00 

0.360 

1.2 

1.1 

14.2 

0.29 

15.00 

15.00 

1.00 

0.165 

0.068 

1.42E  +  01 

2.03E  +  03 

32.73 

79.76 

SP28AAAMC 

1.40E  +  400 

0.336 

1.2 

1.1 

14.2 

0.32 

15.00 

15.00 

1.00 

0.135 

0.052 

1.32E  +  01 

1.89E  +  03 

37.33 

96.98 

Electrochemi  stry 

Part  number 

Discharge 
current  (A) 

Energy  (Wh) 

Nominal 
voltage  (V) 

X 

y 

z 

n 

s 

w 

M  (kg) 

Area  (cm2)  Number  of 
cycles  (not 
including 
recharge 
cycles) 

Lifetime  (h) 

Energy  density 
(Whcm-2) 

Bundle  2 

Lithium-free  thin-film  batter 

- 

5.00E  - 

03 

7.2E  -  01 

3.60 

0.41 

4.72 

12.00 

60 

5 

12 

- 

60 

1.5E  +  02 

21176.47 

7.20E  -  01 
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Table  14 

Comparison  of  Approaches  1,  2,  and  3  with  and  without  CMOS  op-amps  and  novel  voltage  generator  technologies  [28-31] 


Results-exclusion  of  op-amp  technology  (inclusion  of  op-amp  technology) 


Approach 

Battery  type(s) 

Number  of  batteries 

Total  mass  (kg) 

Total  volume  (L) 

Lifetime  (h) 

1 

Lithium  polymer 
battery — prismatic 

15(1) 

0.068  (0.014) 

0.058  (0.012) 

2.4E03  (8.13E02) 

2 

Lithium-free  thin-hlm 
prismatic  lithium  ion 

8(4) 

0.032  (0.014) 

0.028  (0.012) 

8.9E03  (8.9E03)  (mW  range),  979  (979) 
(mW  range),  164  (233)  (W  range) 

3 

Lithium-free  thin-hlm 
prismatic  lithium 
polymer 

65  (61) 

0.073  (0.015) 

0.032  (0.006) 

2.1E04  (9.5E02)  (bundle  1),  4.7E03 
(2.1E04)  (bundle  2) 

Commercial  CMOS  op-amp,  and  high-voltage  generator  characteristics 


Manufacturer 

Part  number 

Supply  voltage 
range  Vcc  (V) 

Open  loop 
gain  A  (dB) 

Operational 

temperature 

(C) 

Areal 

footprint 

(cm2) 

Texas  Instruments  [29] 

OPA725 

4-12 

100 

-40  to  125 

0.09 

Texas  Instruments  [29] 

OPA2725 

4-12 

100 

-40  to  125 

0.16 

Seiko  Instruments  [30] 

S-89110A/89120A 

1. 8-5.5 

80 

—40  to  85 

0.04 

University  of  Michigan  [31] 

High-voltage  generator 

Maximum  voltage 
achieved  >200  V 

>300 

— 

0.16 

op-amps)  to  accommodate  for  voltage  gains  above  nominal 
voltages  of  batteries  in  the  database.  Instead,  system  volt¬ 
ages  greater  than  nominal  voltages  of  the  database  batteries, 
were  accommodated  by  placing  batteries  in  series.  In  this 
section  of  the  document,  we  briefly  examine  implementation 
of  op-amp  or  voltage  generator  technology  to  achieve  volt¬ 
ages  substantially  higher  than  the  nominal  battery  voltage 
values.  We  explore  use  of  op-amps  because  assembling  sep¬ 
arate  transistors,  resistors,  and  capacitors  on  a  circuit  board  is 
an  expensive  and  wasteful  procedure,  while  use  of  integrated- 
circuit  (IC)  amplifiers  is  preferential  because  they  are  versa¬ 
tile,  compact,  and  inexpensive.  Op-amps  are  a  collection  of 
circuit  elements  used  to  increase  voltage,  current,  or  power 
within  a  system,  that  can  be  single  stage  (consisting  of  a  single 
transistor)  or  multistage  (more  than  one  stage);  e.g.  voltage 
follower,  noninverting,  inverting,  summing,  and  current- to- 
voltage  converter.  We  have  selected  noninverting  op-amps 
[29,30]  or  voltage  generators  [3 1],  to  achieve  voltage  gain  for 
the  three  approaches  providing  voltage  ratio  values  (y)  greater 
than  1.  Table  14  details  some  typical  characteristics  of  com¬ 
mercial  CMOS  op-amp  devices  along  with  a  novel  voltage 
generator  technology  demonstrated  by  Udeshi  [31].  Table  14 
summarizes  results  for  all  three  approaches  for  cases  where 
voltage  generators  [31]  are  used  and  not  used.  The  mass  and 
volume  of  the  voltage  generator(s)  are  not  included  in  the 
mass  and  volume  calculations. 


6.  Discussion 

Table  14  summarizes  the  results  from  each  approach.  The 
advantages  of  using  these  approaches  are  that  it  is  a  systematic 
method  for  the  assembly  of  voltage  and  discharge  current  data 
that  readily  lends  itself  to  the  application  of  conditionality 


statements  and  step-wise  selection  techniques  that  deliver  an 
optimal  solution  to  a  general  MEMS  power  problem. 

All  of  these  approaches  produce  mass  and  volume  values 
that  surpass  the  target  values,  0.0043  kg  and  0.00145  cm3. 
The  specific  energy  and  energy  density  values  calculated  for 
each  approach  assume  theoretical  capacity  of  active  materi¬ 
als.  In  fact,  the  actual  energy  available  from  a  real  battery 
under  practical,  but  close  to  optimal  discharge  conditions  is 
only  about  25-35%  of  the  theoretical  energy  of  the  active 
materials  [17],  because  battery  manufactures  must  include 
electrolyte,  packaging,  and  additives  to  enhance  conductiv¬ 
ity  in  the  calculation  of  the  specific  energy  and  energy  density 
of  their  products. 

Approach  2  provided  the  best  results  of  mass  (0.032  kg) 
and  volume  (0.028  L)  among  the  three  approaches.  Specifi¬ 
cally,  use  of  Approach  2  results  in  masses  47.1  and  43.8% 
and  volumes  48.3  and  87.5%  of  Approaches  1  and  3,  re¬ 
spectively.  In  general,  Approaches  2  and  3  provided  the  best 
battery  lifetime  results;  both  systems  produced  lifetimes  in 
excess  of  2E3  h.  Approach  3  required  the  most  batteries  (4.3 
and  8.1  times  Approaches  1  and  2),  which  increases  the  prob¬ 
ability  of  system  failure  and  fabrication  complexity. 

Approach  1  is  appropriate  for  systems  having  relatively 
low  and  constant  current  drain  rates  and  nominal  voltage  val¬ 
ues,  small  to  moderate  cycle  times,  and  sufficient  time  for 
slow  to  moderate  (4-5  h)  charging  time.  Also  this  approach 
provides  additional  benefits  of  a  less  design  and  fabrication 
complexity  since  only  one  electrochemistry  is  used.  How¬ 
ever,  additional  capacity  or  battery  lifetime  achieved  from 
the  additional  cells  for  accommodation  of  fluctuations  in  dis¬ 
charge  current  and  voltage;  can  potentially  result  in  heavy  and 
large  systems.  Extensive  power  management  systems  could 
be  necessary  to  adequately  manage  the  different  discharge 
current  and  voltage  requirements  of  the  system. 
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Approach  2  capitalizes  on  the  variability  of  battery  per¬ 
formance  based  on  shape  and  fabrication,  to  reduce  mass  and 
volume  power  requirements.  It,  unsurprisingly,  resulted  in 
the  smallest  mass  and  volume  values  (0.032  kg  and  0.028  L, 
respectively)  for  the  systems  studied.  Also,  this  strategy  al¬ 
lows  power  size,  shape,  and  manufacturability  to  be  consid¬ 
ered  for  optimal  battery  lifetime  and  system  mass  and  volume 
constraints,  as  demonstrated  with  the  implementation  of  thin- 
film  battery  technology  for  the  micro-power  range,  and  thin 
prismatic  cells  for  both  the  mW  and  W  power  ranges. 

Despite  the  many  advantages  of  Approach  2,  it  is  not  suit¬ 
able  for  all  systems,  for  example,  power  system  fabrication 
complexity,  could  be  increased  as  each  battery  type  could 
require  a  different  fabrication  technique.  This  approach  also 
generally  does  not  account  for  battery  lifetime;  in  the  case 
studied  here,  it  provided  the  shortest  battery  lifetime  pre¬ 
dictions  for  the  mW  and  W  power  ranges  (979  and  164  h, 
respectively,  excluding  the  jjiW  power  range).  Also,  use  of 
different  battery  electrochemistries  requires  different  charg¬ 
ing  programming  for  each  electrochemistry  used. 

Approach  3  provided  the  worse  values  for  mass  (0.073  kg) 
of  all  approaches,  however,  provided  a  volume  very  close  to 
Approach  2  ( — '1.14  times).  However,  use  of  smaller  stand¬ 
alone  power  systems  is  a  safe  guard  against  complete  system 
failure  that  can  occur  when  only  one  power  system  is  used. 
This  approach  allows  an  additional  degree  of  freedom,  in  that 
instead  of  a  maximum  of  only  three  electrochemistries  being 
an  option  for  the  power  system,  the  maximum  number  of 
electrochemistries  that  can  be  implemented  is  equal  to  the 
number  of  subdevices  and/or  power  bundles. 

It  is  obvious  from  comparison  of  results  listed  in  Table  14 
that  use  of  voltage  generator  or  noninverting  op-amp  tech¬ 
nology  can  certainly  reduce  the  mass  and  volume  of  a  power 
system  in  cases  where  power  ranges  are  low  to  moderate, 
while  voltage  values  are  in  excess  of  3.7  V.  However,  mass, 
volume  and  surface  area  of  op-amps  and/or  voltage  genera¬ 
tors  must  also  be  weighed  in  power  system  design,  along  with 
required  power  and  discharge  current  for  operation.  Also,  op- 
amp  performance  as  a  function  of  temperature  is  important  as 
drift  (inaccurate  voltage  gain  values  due  to  temperature  vari¬ 
ation)  becomes  more  probable  at  high-operation  temperature 
temperatures. 

Use  of  multiple  battery  power  sources  can  increase  the 
probability  of  system  failure  due  to  the  increase  in  the  num¬ 
ber  of  fabrication  steps  and  power  subdevices.  In  this  case 
study,  65  batteries  are  required;  other  workers  have  also  im¬ 
plemented  multiple  thin-film  batteries  in  series  and  parallel 
[10].  Also,  this  approach  is  limited  by  additional  mass,  vol¬ 
ume,  and  areal  footprint  that  may  be  necessary  for  individual 
power  supplies  for  each  component  along  with  additional 
manufacturing  time. 

In  this  case  study,  user  prescribed  power  bundle  configu¬ 
rations  were  used  based  on  packaging  and  continuity  of  fluid 
flow  constraints,  which  precluded  detailed  consideration  of 
locations  of  power  bundle  sites.  This  has  a  profound  impact 
on  results  obtained  here,  as  the  WIMS-EMT  consists  of  10 


subdevices,  which  lends  itself  to  numerous  permutations  that 
were  not  examined  here. 

7.  Conclusions 

Through  this  analysis  of  the  WIMS-EMT,  we  have 
demonstrated  that  the  solution  to  power  supply  problem  will 
most  likely  be  a  nonlinear  optimization  problem,  subject  to 
both  linear  and  nonlinear  constraints  and  that  the  most  effi¬ 
cient  power  design  will  most  likely  be  a  hybrid  power  system 
based  on  the  superior  performance  of  both  hybrid  systems 
from  Approaches  2  and  3.  Approach  2  provided  the  best  re¬ 
sults  of  mass  (0.032  kg)  and  volume  (0.028  L)  among  the 
three  strategies.  Specifically,  use  of  Approach  2  results  in 
masses  47.1  and  43.8%  and  volumes  48.3  and  87.5%  of  Ap¬ 
proaches  1  and  3,  respectively.  In  general,  Approaches  1  and 
3  provided  the  best  battery  lifetime  results  wherein  both  pre¬ 
dicted  lifetimes  in  excess  of  2E3  h.  Approach  3  required  the 
most  batteries  (4.3  and  8.1  times  Approaches  1  and  2),  which 
increases  the  probability  of  system  failure  and  fabrication 
complexity. 

Diversity  of  battery  types  within  the  database,  i.e.  high- 
energy  density  and  high-power  density  materials  is  impera¬ 
tive  for  effective  power  supply  design.  The  availability  of  bat¬ 
tery  performance  data,  such  as  cell  capacity  as  a  function  of 
discharge  rate  is  crucial,  as  calculation  of  energy  and  current 
factors  (x  and  y)  depends  on  it.  Specifically,  relationships  be¬ 
tween  discharge  rate  and  capacity  can  prevent  the  elimination 
of  batteries  unnecessarily,  such  as  when  only  one  capacity  at 
one  discharge  rate  is  provided.  Results  from  the  implementa¬ 
tion  of  this  algorithm  are  only  as  reliable  as  the  source  of  the 
battery  characteristics:  the  battery  database.  Incomplete  data, 
i.e.  little  to  no  data  on  capacity  of  battery  at  various  discharge 
rates,  can  result  elimination  of  viable  power  solutions  early 
in  the  optimization  process. 

8.  Future  work 

Our  method  of  predicting  battery  life  was  clearly  approx¬ 
imate.  By  assuming  all  cycles  occur  continuously,  battery  re¬ 
covery,  i.e.  recovery  of  battery  voltage  that  occurs  when  a  bat¬ 
tery  stands  idle  after  a  discharge,  was  not  taken  into  account. 
Also,  potential  battery  self-discharge  was  not  taken  into  ac¬ 
count.  Reduction  in  battery  capacity  with  subsequent  dis¬ 
charge  cycles  was  similarly  not  accounted  for  in  the  present 
work.  Collectively,  these  assumptions  could  result  in  signifi¬ 
cant  overestimates  of  battery  lifetime,  and  will  be  investigated 
as  part  of  future  work. 

Experimental  verification  of  code  predictions  will  also  fol¬ 
low.  Other  possible  selection  criteria  not  investigated  in  this 
document,  but  of  importance,  are  power  system  cost,  number 
of  manufacturing  steps,  complexity,  probability  of  failure  due 
to  system  complexity,  and  temperature.  Of  these,  we  will  next 
focus  on  the  probability  of  failure  due  to  system  complexity, 
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along  with  incorporation  of  some  CMOS  op-amp  technolo¬ 
gies  to  accommodate  large  voltage  fluxes  in  many  MEMS  de¬ 
vices.  Also,  other  means  of  grouping  based  on  power  bundle 
location,  wherein,  varying  the  number  of  devices  within  a  set 
number  of  bundle  locations  will  be  investigated.  Lastly,  work 
in  establishing  an  effective  method  for  prediction  of  battery 
capacity  as  a  function  of  electrochemical  reaction  rates  will 
be  pursued.  Ultimately,  we  aim,  as  part  of  future  work,  to  im¬ 
plement  our  verified  approaches,  conditionality  statements, 
and  selection  techniques  within  a  user-friendly  code. 
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